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Abstract: Ongoing deep surveys of galaxy luminosity distribution functions, spectral energy distribu- 
tions and backwards evolution models of star formation rates can be used to calculate the past history of 
intergalactic photon densities and, from them, the present and past optical depth of the universe to 7-rays 
from pair production interactions with these photons. Stecker, Malkan & Scully have recently calculated 
the densities of intergalactic background light (IBL) photons of energies from 0.03 eV to the Lyman limit 
at 13.6 eV and for redshifts 0< z < 6, using deep survey galaxy observations from the Spitzer; Hub- 
ble and GALEX space telescopes. From these results, they have predicted absorption features for blazar 
spectra. 

This proceedure can also be reversed by looking for sharp cutoffs in the spectra of extragalactic 7-ray 
sources at high redshifts in the multi-GeV energy range with GLAST (the Gamma-ray Large Area Space 
Telescope). Determining the cutoff energies of sources with known redshifts and little intrinsic absorption 
may enable a more precise determination of the IBL photon densities in the past, i.e., the "archaeo-IBL", 
and therefore allow a better measure of the past history of the total star formation rate, including that from 
galaxies too faint to be observed. 

Conversely, observations of sharp high energy cutoffs in the 7-ray spectra of blazars at unknown redshifts 
can be used instead of spectral lines to give a measure of their redshifts. Also, given a knowledge of 
the archaeo-IBL, one can derive the intrinsic 7-ray spectra and luminosities of blazars over a range of 
redshifts and look for possible trends in blazar evolution. Stecker, Baring & Summerlin have found some 
evidence hinting that TeV blazars with flatter spectra have higher intrinsic TeV 7-ray luminosities and 
indicating that there may be a correlation of flatness and luminosity with redshift. GLAST will observe 
and investigate many blazars in the GeV energy range and will therefore provide much new information 
regarding this possibility. 



Introduction 

Space-based telescopes that are dedicated to ex- 
ploring the Universe out to large distances in wave- 
length ranges from the far-infrared to the X-ray 
range are now in place and the Gamma-ray Large 
Area Space Telescope, GLAST, is scheduled to be 
early next year. These facilities are capable of 
probing the Universe to study the early primor- 
dial phases of star formation and galaxy evolution. 
Studies of cosmic photons from NASA space ob- 
servatories supply important inputs as well into as- 
trophysics, physics and cosmology. Presently, the 
Spitzer space infrared telescope facility is probing 
deeply into the past to study galaxy formation and 
evolution in the infrared. GALEX is making sim- 
ilar observational studies in the ultraviolet. In the 



7-ray energy range SWIFT is gathering important 
information on 7-ray bursts back to the distant past 
and GLAST will study sources and diffuse fluxes of 
7-rays at energies up to 100 GeV. The astronomical 
era of the "deep fields" has arrived. 

Whereas archaeological sites such as Chi-chen 
Itza, which many of us visited on the conference 
excursion, allow us to peer back to '^kyr in the 
past, the deep astronomical surveys from the ra- 
dio to the 7-ray wavelengths allow us to peer back 
^Gyr or more into the history of the Universe. 



Photon and Blazar Archaeology 





lOOum 




10um 


40 


, 1 1 1 11 1 1 1 1 


, 1 , , 


, 1 1 1 1 . 
1000L* 


38 


/ 




- 


g 


/ 
















->-. 






s 

3 








J* 36 






34 


, 1 . . , , 1 . . . 


, 1 , , 


.001 L* 



12 12.5 13 13.5 14 

Log Frequency (Hz) 



Figure 1 : Empirically derived galaxy SEDs for var- 
ious luminosities from 10~^L* to lO^i*. 

Calculating Intergalactic Archaeopho- 
ton Densities 

Stacker, MaUcan & Scully [1] have used the ap- 
proach pioneered in Refs. [2] and [3] to calculate 
past intergalactic infrared photon fluxes and den- 
sities as they evolved along with the galaxies that 
produced them. This method, a "backwards evo- 
lution" scheme, is an empirically based calculation 
which uses as input (1) the luminosity dependent 
galaxy spectral energy distributions (SEDs) based 
on observations of normal galaxy infrared SEDs, 
(2) observationally based galaxy luminosity distri- 
bution functions (LPs) and (3) redshift dependent 
galaxy luminosity evolution functions. These are 
empirically derived curves giving the universal star 
formation rate or luminosity density as a function 
of redshift, which are referred to as Lilly-Madau 
plots. Ref. [1] uses new deep survey observations, 
primarily from the Spitzer infrared space observa- 
tory, to make improvements to the previous calcu- 
lations of Refs. [2] and [3]. 

Galaxy Infrared SEDs as a Function of Lu- 
minosity 



mines the average galaxy SED and that therefore 
the galaxy luminosity distribution function (LP) 
can be predicted statistically from its observed lu- 
minosity in one infrared waveband, here chosen 
to be 60/im [4J. It is a well-established fact that 
more luminous galaxies (now and in the past) have 
higher rates of ongoing star formation and that the 
star formation rate was higher in the past. 

Empirically, it is also found that for the more lumi- 
nous galaxies, relatively more of the energy from 
these young stars is absorbed by dust grains and 
re-radiated in the thermal infrared; more luminous 
galaxies have higher infrared flux relative to opti- 
cal flux and warmer infrared spectra. These clear 
luminosity-dependent trends in galaxy SEDs were 
well determined locally from the combination of 
IRAS (Infrared Astronomy Satelhte) and ground- 
based photometry for large (e.g., all sky) samples. 
The infrared SEDs as a function of galaxy luminos- 
ity used in [1] were based on broadband photome- 
try of IRAS selected samples and show the average 
SED emission trends discussed above. The fam- 
ily of average galaxy SEDs at various luminosities 
in terms of L* is shown in Figure 1, taken from 
Ref. [2], where = 8.4 x 10^^ W Hz'^ is flie 
luminosity at 60 /xm. 

Other computations of infrared backgrounds and 
source counts have used different SEDs, based on 
somehwat different combinations of data and mod- 
els, in some cases estimated in more spectral detail. 
One may ask whether these new SEDs might dif- 
fer from those used in Refs. [2], [3] and [I], either 
in overall colors, or in detail around the 7-12/Ltm 
region, where the strongest spectral features from 
polycyclic aromatic hydrocarbon (PAH) emission 
and siUcate emission are found. The best example 
of these new galaxy template SEDs can be found 
in Ref [5]. As shown in Ref. [I], the agreement 
between the infrared SEDs in Refs. [5] and [1] 
is excellent. As long as this agreement holds for 
SEDs of galaxies near the "knee" of the approxi- 
mate broken power-law LFs is reasonable (see Fig- 
ure 2), the final computed intergalactic photon den- 
sities will also agree, because they are dominated 
by galaxies with luminosities ~ i# around the 
knee. 



The key empirically supported assumption made in 
Ref. [1] is that the luminosity of a galaxy deter- 
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Figure 2: 60/im local galaxy luminosity function. 

The Local Infrared Luminosity Function 

The foundation of the backwards evolution calcu- 
lation [ 1 ] is an accurately determined local infrared 
luminosity function of galaxies. We used the LF 
observationally derived in Ref . [4] because it was 
based on an extensive analysis of very large data 
sample of galaxies. This LF was updated by using 
the even more thorough local infrared LF given in 
the more recent Ref. [6]. A low-luminosity power- 
law index of -1.35 was adopted in Ref. [1] for the 
broken power-law differential LF which steepened 
by 2.25 at high luminosities. This LF takes bet- 
ter account of the large number of fainter galaxies 
that are now known [7]. The cosmological values 
adopted in Ref. [1] were h = 0.7 and a ACDM 
cosmology with = 0.7 and ^Im = 0.3. 

Figure 2 compares the broken power-law LF at 
z = with the data given in Ref. [6]. Note that 
the luminosity used in Figure 1 is defined to be 
the luminosity at the knee of the LF shown in Fig- 
ure 2. It can be seen that the agreement between 
the analytic curve and the observational data is ex- 
cellent. 

Evolution of the IBL SED with Redshift 

It is now well known that galaxies had a brighter 
past owing to their higher rates of star formation 



and the fading of stellar populations as they age. 
The simplest resulting evolution of the galaxy lu- 
minosity function is a uniform shift in either the 
vertical axis (number density evolution), or in the 
horizontal axis (luminosity evolution.) For a pure 
power-law luminosity function, number and lumi- 
nosity evolution are mathematically equivalent. In 
reaUty, however, to avoid unphysical divergences 
in the total number or luminosity of galaxies, the 
luminosity function must steepen at high L and 
flatten at low L. Thus real LF's will have at least 
one characteristic "knee" separating the steep high- 
L portion from the flatter low-L slope. For typical 
LF's this results in most of the luminosity being 
emitted by galaxies within an order of magnitude 
of this knee. Thus large uncertainties and errors in 
the LF far from this knee wiU hardly change most 
of the results (e.g., number counts and integrated 
diffuse backgrounds). 

Strong luminosity evolution of galaxies, i.e., a sub- 
stantial increase in the luminosity of this knee 
with redshift, is consistently found by many ob- 
servations relating infrared luminosity to the much 
higher star formation rate at 2; ~ 1 and to the recent 
determination that most ultraviolet-selected galax- 
ies at z ~ 1 are also luminous infrared galaxies. 

The exact form of the luminosity evolution of 
galaxies is the most uncertain input to the calcula- 
tion of the IBL density, since it depends on high 
redshift data. Therefore, two plausible cases of 
pure luminosity evolution were adopted in Ref. [1] 
in order to bracket the uncertainties from the deep 
surveys: (1) a more conservative "baseline" (B) 
evolution scenario and (2) a "fast evolution" (FE) 
scenario. A more detailed discussion of these mod- 
els is given in Ref. [1]. The FE model is favored 
by recent Spitzer observations [8, 9]. It also pro- 
vides a better description of the deep Spitzer num- 
ber counts at 70 and 160/im than the B model. 
However, GALEX (Galaxy Evolution Explorer) ob- 
servations indicate that the redshift evolution of ul- 
traviolet radiation may be more consistent with the 
B model [10]. The Spitzer IRAC (Infrared Array 
Camera) counts can be best fit with an evolution 
rate between these two models. One way of under- 
standing the somewhat smaller redshift evolution 
of the star formation rate implied by the GALEX 
ultraviolet observations vs. that obtained from the 
Spitzer infrared observations is that the effect of 



Photon and Blazar Archaeology 



dust extinction followed by infrared reradiation in- 
creases with redshift [11 J. 



Calculation of the Intergalactic Back- 
ground Light 

Our calculation of the diffuse infrared background 
as described above extends up to a rest frequency 
of log uhz = 14.1. corresponding to an energy 
of ~ 0.5 eV. This is the location of the peak in 
the spectral energy distributions of most galax- 
ies and is produced by the light of red giant 
stars. The spectrum (in energy density units) then 
curves downward rapidly to higher frequencies, 
with modest dependence on galaxy luminosity. Al- 
though galaxy SEDs have a peak at this energy, this 
peak, as opposed to the far infrared peak in galaxy 
SEDs, only manifests itself as an inflection point 
in the photon density spectrum.^ At wavelengths 
shortward of this near infrared-optical "peak' the 
photon density spectra drop steeply (See Figures 
3 and 4 which plot en(e) for the case of the fast 
evolution (FE) model.) The relative increase in the 
ratio of red-to-blue photons at the higher redshifts 
shown in these figures is produced by early gener- 
ation, low luminosity stars. 

The low redshift extragalactic background light 
(EBL) SEDs obtained are shown in Figure 5. They 
are given for the two evolution models previously 
discussed, viz., the "baseline" (B) model and the 
"fast evolution" (FE) model. 



Gamma-ray Absorption from Pair Pro- 
duction Interactions 
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Figure 3: The photon density en{e) shown as a 
continuous function of energy and redshift. 
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The potential importance of the photon-photon pair 
production process, 77 e+e~, in high energy 
astrophysics has been realized for over 40 years 
[14]. It was pointed out that owing to interactions 
with the 2.7 K CMB, the universe would be opaque 
to 7-rays of energy above 100 TeV at extragalac- 
tic distances [15, 16]. If one considers cosmo- 
logical and redshift effects, it was further shown 
that photons from a 7-ray source at a redshift Zs 
would be significantly absorbed by pair produc- 
tion interactions with the CMB above an energy 
~ 100(1 + 2^)"^ TeV [17, 18]. 



Figure 4: The photon density en(e) function of en- 
ergy for various redshifts based on the fast evolu- 
tion model for infrared evolution. The solid line is 
for z = 0, the dashed line is for 2; = 1 , the dotted 
Une is for 2; = 3, the dot-dashed Une is for z = 5. 



1 . We note that the energy dependence of the differ- 
ential photon energy spectrum, dn~,/de, is obtained by 
dividing the SED by the square of the photon energy 
(e^) so that the starlight "peak" in the SED has very few 
photons compared to the dust reradiation peak in the far 
infrared. 
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Figure 5: Spectral energy distribution of the diffuse 
background radiation at 2 = 0, sometimes referred 
to as the extragalactic background light (EBL). Er- 
ror bars show data points, triangles show lower 
limits from number counts, and the inverted trian- 
gle shows an upper limit from TeV observations 
[12]. The upper and lower sohd lines show the 
fast evolution and baseline evolution predictions, 
and the dotted lines show their extensions into the 
optical-ultraviolet, as described in Ref. [13]. 



Following the discovery by the EGRET team of the 
strongly flaring 7-ray blazar 3C279 at redshift 0.54 
[19], Stecker, de Jager and Salamon [20] proposed 
that one can use the predicted pair production ab- 
sorption features in blazars to determine the inten- 
sity of the infrared portion of the IBL, provided 
that the intrinsic spectra of blazars extends to TeV 
energies. It was later shown that the IBL produced 
by stars in galaxies at redshifts out to 2 would 
make the universe opaque to photons above an en- 
ergy of ~ 30 GeV emitted by sources at a redshift 
of ^ 2, again owing to pair production interactions 
[21, 13]. 

As discussed above, in Ref. [1] this approach was 
expanded by using recent data from the Spitzer in- 
frared observatory In addition, data from the Hub- 
ble deep survey and the GALEX mission were also 
used to determine the photon density of the IBL 
from 0.03 eV to the Lyman limit at 13.6 eV for 
redshifts out to 6 (the "archaeo-IBL").^ The results, 
giving the IBL photon density as a function of red- 
shift together with the opacity of the CMB as a 
function of redshift, were then used to calculate the 
opacity of the universe to 7-rays for energies from 
4 GeV to 100 TeV and for redshifts from - to 5. 

The results of Ref. [1] as shown in Figure 6 imply 
that the universe will become opaque to 7-rays for 
sources at the higher redshifts at somewhat lower 
7-ray energies than those given in Ref. [13]. This 
dependence is shown in Figure 8. The increased 
opacity calculated at the higher redshifts is because 
the newer deep surveys have shown that there is 
significant star formation out to redshifts z > 6 
[25, 26]), greater than the value of z^ax = 4 as- 
sumed in Ref. [13]. This conclusion is also sup- 
ported by recent Swift observations of the redshift 
distribution ofGRBs [27]. 

Stecker, Malkan and Scully [1] found that the func- 
tion t{E^,z) shown in Figure 6 can be very well 
approximated by the analytic form 



logT = Ax^ + Bx^ + Cx"^ + Dx + E (1) 

over the range 0.01 < r < 100 where x = log 
(eV). The correct coefficients A through E are 
given for various redshifts in Ref . [28]. (This is the 



2. See also Refs. [22,23,24]. 
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Table 1. Optical Depth Parameters for the Fast Evolution and Baseline Models 



Evolution Model 


A 


B 


C 


D 


Fast Evolution 


-0.475 


21.6 


-0.0972 


10.6 


Baseline 


-0.346 


16.3 


-0.0675 


7.99 




le+OO \e*0\ le+Q: lc+U3 le+tl4 le+05 



Figure 6: Optical depth of the universe to 7-rays 
from interactions with photons from the EBL and 
2.7K CMB for various redshifts, z. The solid lines 
are for the EE model and the dashed lines are for 
the B model. The curves shown are for (from top 
to bottom) z = 5.0, 3.0, 2.0, 1.0, 0.5, 0.2, 0.117, 
0.03[1]. 



corrected form of the original table given in Ref. 
[1].) 
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Figure 7: Schematic of the GLAST satellite de- 
ployed in orbit. The LAT, which comprises of a 
cluster of 16 silicon strip trackers and a calorime- 
ter consisting of 8 layers of CsI(Tl) crystals, is in 
the top (yellow) area; the Gamma Ray Bust Moni- 
tors GBM are located directly below. 

Observations of relatively nearby blazars using the 
H.E.S.S. and MAGIC air Cerenkov telescopes have 
produced many results in the TeV energy range. 
However, as can be seen from Figure 6, the 
critical energy range for exploring sharp absorp- 
tion cutoffs from distant blazars is the multi-GeV 
range. This is an energy range which the up- 
coming GLAST mission will explore. The GLAST 
satellite (see Figure 7) will have a large area 
telescope called the LAT which is designed to 
study cosmic diffuse 7-rays and 7-ray sources in 
the energy range between '^20 MeV and > 300 
GeV (For more information about GLAST, go to 
http://glast.stanford.edu/. By using 
GLAST observations to determine the absorption 
cutoff energies of sources with known redshifts 
caused by interactions of GeV range 7-rays with 
low energy photons of the IBL, we can refine our 
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knowledge of the IBL photon densities in the past, 
i.e., the archaeo-IBL, and therefore get a better 
measure of the past history of the total star forma- 
tion rate. Conversely, observations of sharp high 
energy cutoffs in the 7-ray spectra of sources at 
unknown redshifts can be used instead of spectral 
lines to give a measure of their redshifts. GLAST 
must also investigate the possibility that intrinsic 
absorption within some high luminosity, high red- 
shift 7-ray sources may mimic absorption by the 
IBL [29]. If this is the case, such an effect will 
need to be disentangled from IBL absorption. One 
key question is how far the region of 7-ray emis- 
sion is from the radiation field surrounding the 
black hole. There evidence in the case of M87 that 
the 7-ray emission region is greater than 120 pc 
from the black hole [30] arguing against intrinsic 
absorption. However, this is only one source and 
it is not a typical blazar. Observationally, intrin- 
sic absorption in active galactic nuclei can be in- 
vestigated by determining the lower limit envelope 
on the total absorption from a set of sources in a 
known narrow redshift range. This can be accom- 
pUshed by GLAST. 

Blazar Spectra and Luminosities, Now 
and Then 

As an example of the application of the 7-ray ab- 
sorption models derived in Ref. [1], Figure 9 shows 
the absorbed spectrum of the blazar PKS 2155-304 
at z = 0.117 assuming that its intrinsic spectrum 
in the energy range shown is an approximate E~'^ 
power law. This spectrum is then compared it with 
the spectrum observed by the H.E.S.S. TeV 7-ray 
telescope array. Thus, one would conclude that 
the intrinsic (unabsorbed) photon spectrum of this 
source in the energy range shown has an approxi- 
mate form. 

The intergalactic 7-ray absorption coefficient (i.e. 
optical depth), t{E, z), increases monotonically 
with energy and therefore leads to a steepening 
of the intrinsic source spectra as it is observed 
at Earth. For sources at redshifts between 0.05 
and 0.4, Stecker & Scully [3 1] have shown that 
this steepening results in a well-defined increase 
in the spectral index of a source with an approxi- 
mate power-law spectrum in the 0.2 - 2 TeV energy 
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Figure 8: A Fazio-Stecker plot (named after the 
paper by Fazio & Stecker [18]) which gives the 
critical optical depth r = 1 as a function of 7- 
ray energy and redshift for the fast evolution (soUd 
curve) and baseline (dashed curve) IGL cases. Ar- 
eas to the right and above these curves correspond 
to the region where the universe is optically thick 
to 7-rays (from Ref. [1].) 
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Figure 9: The 7-ray data from the H.E.S.S. for 
PKS2155-304 compared with the theoretically ob- 
tained spectrum for PKS2155-304 calculated by 
assuming an intrinsic source spectrum proportional 
to and multiplying by using t{z = 
0.117) for the FE model of Ref. [1]. 
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range. This increase is a linear function in redshift 
z of the form Ar = C + Dz, where the parame- 
ters C and D are constants. The overall normaliza- 
tion of the source spectrum is also reduced by an 
amount equal to exp{ — {A + Bz)} , again where 
A and B are constants. The values of A, B, C, and 
D are given for the B and FE models in Table 1 . 

These analytic relations can be to calculate the in- 
trinsic 0.2 - 2 TeV power-law 7-ray spectra of 
sources having known redshifts in the 0.05 - 0.4 
range and observed spectral indeces, Fq for both 
the B and FE models of EBL evolution Table 2, 
taken from Ref. [32], gives values for Tg for var- 
ious blazars obtained using the formula given in 
Ref. [31]. Table 1 also shows the respective in- 
dices Fe = 2Fs — 1 of the electron distributions 
in the sources under the assumption that the 7-rays 
are produced by inverse Compton interactions in 
the Thomson regime. 

Stecker, Baring & Summerlin [32] have also esti- 
mated the intrinsic "isotropic luminosity." of these 
blazars, also shown in Table 2^ The isotropic lu- 
minosity of the blazar sources listed in Table 2 is 
obtained from the formula 

jC ^ 4n^^{l+zf'-'FMz)fe^^^'''^ (2) 

where d is the luminosity distance to the source, 
and Fo is its observed differential energy flux at 

1 TeV, and the other factors in the equation give 
the k-correction for the deabsorbed source spec- 
trum and the normalization correction factor for 
absorption, exp{ — {A + Bz)} given in Ref. [31 J. 
The blazars Mrk 421 and Mrk 501 are not in- 
cluded because their redshifts are significantly less 
than 0.05. However, these blazars are analysed in 
Ref. [33]. 

The numbers given in the last column of Table 

2 are derived for the fast evolution (FE) model. 
One may note that there appears to be a trend 
toward blazars having flatter intrinsic TeV spec- 
tra and higher isotropic luminosities at higher red- 
shifts. However, one must be careful of selection 
effects. The TeV photon fluxes of these sources as 
observed by H.E.S.S. and MAGIC only cover a dy- 
namic range of a factor of ~20. Therefore, only 
brighter sources can be observed at higher red- 
shifts. This is because of both diminution of flux 
with distance and intergalactic absorption [20]. 



The observed luminosity-redshift trend is naturally 
expected in a limited population sample spanning 
a range of redshifts if the TeV-band fluxes are 
pegged near an instrumental sensitivity threshold. 
A more powerful handle on the intrinsic spectra 
and luminosities of these sources will be afforded 
by the upcoming GLAST 7-ray mission, with its 
capability for detecting many blazars at energies 
below 200 GeV. 

Table 2 indicates that some blazars, particularly 
those observable at the higher redshifts, appear to 
have very hard intrinsic 7-ray spectra with indeces 
in the range between ^^1 and '--^1.5. The possibil- 
ity that such spectra can be obtained from shock 
acceleration has not usually been admitted when 
considering properties of blazar jets and their pos- 
sible emission spectra [34]. However, such spec- 
tra can be produced by relativistic shock accelera- 
tion [32], as discussed in the next section. We will 
also discuss observational evidence in the hard X- 
ray range for 7-ray blazars with spectral indeces in 
the range between ~1 and ~1.5. 

Shock Acceleration and Intrinsic 
Blazar Spectra 

The rapid variability seen in TeV flares drives the 
prevailing picture for the blazar source environ- 
ment, one of a compact, relativistic jet that is struc- 
tured on small spatial scales that are unresolvable 
by present 7-ray telescopes. Turbulence in the su- 
personic outflow in these jets naturally generates 
relativistic shocks, and these form the principal 
sites for acceleration of electrons and ions to the 
ultrarelativistic energies imphed by the TeV 7-ray 
observations. Within the context of this relativistic, 
diffusive shock acceleration mechanism, numeri- 
cal simulations can be used to derive expectations 
for the energy distributions of particles accelerated 
in blazar jets. 

Diffusive acceleration at relativistic shocks is less 
well studied than for nonrelativistic flows, yet it 

3. We define isotropic here as if the source had an 
apparent isotropic luminosity even though blazars are 
highly beamed and their flux (and hence their appar- 
ent luminosity) is dramatically enhanced by relativistic 
Doppler boosting. This is similar to the nomenclature 
used for 7-ray bursts. The quantity C is equal to AirvFi, 
given at hi^ = 1 TeV in units of 10^® W (10*^ erg s'^). 
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Table 2. Blazar Spectral Indeces in the 0.2-2 TeV Energy Range and TeV Luminosities 



Source Name 


Z 


ohs 


r 

^ source 


(FB B) 


r^(FB 




C(l TeV) [10^*^ W 


lES 2344+514 


0.044 


3.0 


2.5 


^2.6 


4.0^ 


4.2 


2.9 


Mrk 180 


0.045 


3.3 


2.9 


3.0 


4.8 ^ 


5.0 


1.2 


IBS 1959+650 


0.047 


2.7 


2.3 


^2.4 


3.6 ^ 


3.8 


5.4 


PKS 2005-489 


0.071 


4.0 


3.4 


^ 3.5 


5.8 


6.0 


8.6 


PKS 2155-304 


0.117 


3.3 


2.2 


^2.4 


3.4 ^ 


3.8 


420 


H 2356-309 


0.165 


3.1 


1.5 


^ 1.9 


2.0- 


2.8 


200 


lES 1218+30 


0.182 


3.0 


1.2 


^ 1.6 


1.4^ 


2.2 


310 


lES 1101-232 


0.186 


2.9 


1.0 


1.5 


1.0^ 


2.0 


230 


IBS 0347-121 


0.188 


3.1 


1.2 


1.7 


1.4^ 


2.4 


1200 


IBS 1101+496 


0.212 


4.0 


1.8 


^2.4 


2.6 ^ 


3.8 


930 



is the most applicable process for extreme objects 
such as pulsar winds, jets in active galactic nu- 
clei, and 7-ray bursts. A key characteristic that 
distinguishes relativistic shocks from their non- 
relativistic counterparts is their inherent anisotropy 
due to rapid convection of particles through and 
away downstream of the shock. This renders an- 
alytic approaches more difficult for ultrarelativis- 
tic upstream flows, though advances can be made 
in special cases, such as the limit of extremely 
small angle scattering (pitch angle diffusion). Ac- 
cordingly, complementary Monte Carlo techniques 
have been employed for relativistic shocks in a 
number of studies, including test-particle analy- 
ses for steady-state shocks of parallel and oblique 
magnetic fields (see refs. in Ref. [32]). This 
approach was used in Ref. [32] to determine key 
spectral characteristics for particles accelerated to 
high energies at relativistic shocks that are of 
relevance to acceleration in blazar jets. For a 
discussion of relativistic shock acceleration, see 
Ref. [35]. 

Representative particle momentum distributions 
that result from the simulation of diffusive accel- 
eration at relativistic shocks are depicted in Fig- 
ure 10 [32]. These distributions are equally appli- 
cable to electrons or ions, and so the mass scale is 
not specified. These results highlight several key 
features. The momentum distributions are flatter 
when produced by faster shocks with a larger up- 
stream flow bulk Lorentz factor of the jet, 71, given 
a fixed velocity compression ratio. This is a con- 
sequence of the increased kinematic energy boost- 
ing occurring at relativistic shocks. Such a charac- 



-o _e 

u 

B 





1 ' 1 


' 1 ' 1 






r = 3, = 0° ~ 






\7A = 30 




























8.C.11 S (LAS) 












9„.„S 30° ^7,^ 






e.c.uS8° _ 






1,1 


1 vWi- , ^ 



2 4 6 8 10 

Logi„[p/(me)] 



Figure 10: Particle distribution functions dN/dp 
from parallel relativistic shocks obtained from a 
Monte Carlo simulation of particle diffusion and 
gyrational transport. Scattering off hydromagnetic 
turbulence is modeled by randomly deflecting par- 
ticle momenta by an angle within a cone, of half- 
angle Bscatt, whose axis coincides with the parti- 
cle momentum prior to scattering. The thick (blue) 
lines are for the large angle scattering cases (LAS). 
These asymptotically approach the power-laws in- 
dicated by thin lines, at high and very high energies 
(not shown). Three smaller angle scattering cases 
are also shown, corresponding to pitch angle dif- 
fusion (PAD). These have high-energy asymptotic 
power-law indices of 1.65, 1.99 and 2.20. 
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teristic is evident for much larger angle scattering, 
as found in the work of other authors referenced 
in Ref. [32]. What is much more striking in Fig- 
ure 10 is that the slope and shape of the nonthermal 
particle distribution depends on the nature of the 
scattering. The asymptotic, ultrarelativistic index 
of Fe = 2.23 is realized only in the mathematical 
hmit of small (pitch) angle diffusion (PAD), where 
the particle momentum is stochastically deflected 
on arbitrarily small angular (and therefore tempo- 
ral) scales. In practice, PAD results when the maxi- 
mum scattering angle 9 scan, is less than the Lorentz 
cone angle I/71 in the upstream region. In such 
cases, particles diffuse in the region upstream of 
the shock only until their angle to the shock normal 
exceeds ~ I/71. Then they are rapidly swept to 
the downstream side of the shock. The energy gain 
per shock crossing cycle is then roughly a factor of 
two, simply derived from relativistic kinematics. 

To contrast these power- law cases. Figure 10 also 
shows the results given in Ref. [32] for large angle 
scattering scenarios (LAS, with 9scait ^ t^), where 
the spectrum is highly structured and much flat- 
ter on average than p~^. The structure, which be- 
comes extremely pronounced for large 71, is kine- 
matic in origin, where large angle deflections lead 
to the distribution of fractional energy gains be- 
tween unity and 712 in successive shock transits by 
particles. Gains like this are kinematically anal- 
ogous to the photon energy boosting by inverse 
Compton scattering. Each structured bump or 
spectral segment shown in Figure 10 corresponds 
to an increment in the number of shock crossings, 
successively from 1 ^ 3 ^ 5 ^ 7 etc. They 
eventually smooth out to asympotically approach 
power-laws that are indicated by the lightweight 
lines in the Figure. The indices of these asymp- 
totic results are all in the range Fg < 2. Inter- 
mediate cases are also depicted in Figure 10, with 
9scatt ^ 4/71. The spectrum is smooth, like the 
PAD case, but the index is lower than 2.23. As- 
trophysically, there is no reason to exclude such 
cases. From the plasma point of view, magnetic 
turbulence could easily be sufficient to effect scat- 
terings on this intermediate angular scale, a con- 
tention that becomes even more salient for ultra- 
relativistic shocks with 71 3> 10. It is also ev- 
ident that a range of spectral indices is produced 
when Oscatt is of the order of I/71. In this case, 



the scattering processes corresponds to a transition 
between the PAD and LAS limits. 

The impUcations of the numerical simulations 
given in Ref. [32] for distributions of relativistic 
particles in blazars are apparent. There can be a 
large range in the spectral indices Fg of the parti- 
cles accelerated in relativistic shocks, and these in- 
dices usually differ from Fg ~ 2.23. They can be 
much steeper, particularly in oblique shocks. How- 
ever, they can also be much flatter, so that quasi- 
power-law particle momentum distributions p^^' 
with Fg < 2 are readily achievable. These spectra 
can be preserved, provided that the electron spectra 
produced during blazar flares are not significantly 
affected by coohng by synchroton radiation. This 
requires that tacc < tccoi which constrains both 
the magnetic field strength and the electron energy 
for gyrorsonant acceleration processes [36]. These 
electrons can then Compton scatter to produce 7- 
ray spectra with F < 1.5. 

In fact, recent obsevations of an extreme MeV 7- 
ray blazar at a redshift of ^3 by Swift |37| and 
the powerful 7-ray quasar PKS 1510-089 at a red- 
shift of 0.361 by Swift and Suzaku [38] both ex- 
hibited power-law spectra in the hard X-ray range 
which were significantly harder than 1.5, implying 
electron spectra significantly harder than value of 2 
usually considered for shock acceleration. Spectra 
in the hard X-ray range do not suffer intergalactic 
absorption so that there is no ambiguity concerning 
their spectral indeces. 

Blazar Archaeology 

It can be seen from Figure 6 that for 7-ray sources 
at the higher redshifts there is a steeper energy de- 
pendence of the optical depth t{E^) near the en- 
ergy where r = 1. There will thus be a sharper 
absorption cutoff for sources at high redshifts. It 
can easily be seen that this effect is caused by the 
sharp drop in the ultraviolet photon density at the 
Lyman limit. 

It is expected that GLAST will be able to resolve 
out thousands of blazars from the general extra- 
galactic 7-ray background [39]. Because of the 
strong energy dependence of absorption in blazar 
spectra at the higher redshifts in the multi-GeV 
range, GLAST wiU be able to probe the archaeo- 
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IBL and thereby probe the early star formation 
rate. GLAST should be able to detect blazars at 
known redshifts z ^ 2 at multi-GeV energies and 
determine their critical cutoff energy. A simple 
observational technique for probing the archaeo- 
IBL has been proposed in Ref. [40]. In such 
ways, GLAST observations at redshifts z > 2 and 
Ej ~ 10 GeV may complement the deep galaxy 
surveys by probing the total star formation rate, 
even that from galaxies too faint to be detected in 
the deep surveys. Future GLAST observations in 
the 5 to 20 GeV energy range may also help to pin 
down the amount of dust extinction in high-redshift 
galaxies by determining the mean density of ultra- 
violet photons at the higher redshifts through their 
absorption effect on the 7-ray spectra of high red- 
shift sources. If the diffuse 7-ray background radi- 
ation is from unresolved blazars [39], a hypothesis 
which can be independently tested by GLAST [41], 
absorption will steepen the spectrum of this radia- 
tion at 7-ray energies above ^ 10 GeV [13]. Thus, 
GLAST can also aquire information about the evo- 
lution of the IBL in this way. 

Conversely, observations of sharp high energy cut- 
offs in the 7-ray spectra of sources at unknown red- 
shifts can be used instead of spectral Unes to give a 
measure of their redshifts. 

Table 2 indicates that there may be an evolution of 
TeV blazars with redshift, both in luminosity and 
spectral index. It may be that past blazars had jets 
with higher bulk Lorentz factors than low-redshift 
BL Lacs. It may also be that these apparent trends 
are only the result of selection effects. There is 
only a small sample of ten sources listed in Ta- 
ble 2. GLAST will be able to explore this ques- 
tion, as it has the potential of seeing thousands of 
blazars [39]. 

Extragalactic Astronomy: The Greatest 
Dig in the Universe 

In keeping with the cosmic high energy theme of 
the 30th International Cosmic Ray Conference, the 
emphasis in this paper has been on high energy 
photons (7-ray s). However, we have seen the in- 
timate connection between these photons and the 
low energy photons that are the subjects of more 
traditional astronomy. In fact, radio surveys, ongo- 



ing infrared-to-ultraviolet deep astronomical sur- 
veys, and surveys of the more newly discovered 
extragalactic 7-ray sources and 7-ray bursts, are 
most effectively used when considered together to 
create a consistent and compelling picture of the 
past history of the Universe. 

Epilogue 

In memory of Frank Culver Jones, my onetime 
mentor and friend for four decades. Frank passed 
away earher this year. He only missed one ICRC 
conference since 1961. He also served as secretary 
to the ICRC. He will be missed by the cosmic-ray 
community. 
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